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By introducing a set of auxiliary equations representing a many-body system, we have derived 
an extension of the Kohn-Sham scheme for the density functional theory. These equations consist 
of a Kohn-Sham-type equation determining single-particle orbitals and an eigen-value equation 
for an effective many-body problem. A variational method similar to the Kohn-Sham technique 
was utilized to derive effective interactions as well as effective potentials without artificial sub- 
stitution of a Hubbard-type interaction and a mean-field correction in the energy functional. 
The second equation is described by an effective many-body Hamiltonian with both 2-body in- 
teractions and mean-field terms. Rigorous formulation of the extended Kohn-Sham equation 
is also given in accordance with the Hadjisavvas-Theophilou formulation. Our formulation can 
be interpreted as a way to define models of the strongly correlated electron systems, e.g. the 
Hubbard model. 
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§1. Introduction 

The density functional theory (DFT) is a general the- 
ory to describe nnn-jjiiiforni interacting electron gas at 
zero tenipeijatiirdll'atra'B^ or in equilibrium at finite tem- 
peratures.oQ'El' This general formalism provides varia- 
tional principle to search the ground state or the density 
matrix of the many-body system through minimization 
of a functional of the single-particle density. It would be 
allowed to say that the modern theory of the electronic 
band-structure calculations would not be established, if 
DFT were not known. DFT is a rigorous theory stating 
that the ground state energy iSjJibtained as a minimum 
of an energy density functional.EfB' EP 

To perform realistic calculation, we have to utilize a 
few steps of reformulation with assumptions and approx- 
imations starting from the rigorous description. The first 
strategy is the Kohn-Sham (KS) scheraej-which provides 
us with an effective one-body problem.oll3) Namely, one 
can obtain an effective equation called the Kohn-Sham 
equation (KSE), which looks like the Schrodinger equa- 
tion for a non-interacting system. Utilizing the so-called 
local-density approximation (LDA)E l£2 L.or the general- 
ized gradient approximation (GGA),^!!^ KSE becomes 
a tractable self-consistent equation with a Hartree term 
and a so-called exchange-correlation potential. 

The DFT-LDA scheme has been regarded as a ma- 
jor strategy of the band-structure calculation. A lot of 
studies of the electronic states have been done to re- 
veal that the KS scheme supplies a well-controlled the- 
oretical description of various material^ including semi- 
conductors, ionic crystals and metals,© as well as the 
large scale systems by virtue of the iterative minimiza- 
tion tedriii(mellj'E£P and description by the localized or- 
bitals .SbS) 



However, it has been recognized that the band calcu- 
lation using DFT-LZMj-does not reproduce several kinds 
of correlation effects,Bc3' e.g. the Mott metal-insulator 
transition and the Kondo effect. Usage of LDA (or GGA) 
has thought to be a serious reason for this difficulty. But, 
the KS scheme itself might not be suitable to describe 
these correlation effects. This is because the density is 
given by an auxiliary non-interacting system. For non- 
interacting Fermions, the ground state is given by a sin- 
gle Slater determinant. This wave function cannot de- 
scribe the correlation effects. Especially, for correlation 
effects originating from short range repulsion, the diffi- 
culty can lead to qualitative errors. 

Since DFT tells us only the single-particle density of 
the ground state, one might say that it is by no means 
hopeful to expect DFT as a general theory to reproduce 
every properties of materials. However, the present KS 
scheme has been practically proven to be a useful skill for 
materials scientists, even when they study single-particle 
excitations, i.e. the band structure. Thus, we intended 
to find an extension of the KS approach keeping its con- 
venience but enlarging a range of applications. 

We expected that a natural extension of the KS scheme 
is to utilize a density given by a many-body state of a 
simplified auxiliary problem. If this scheme becomes pos- 
sible, not only single-particle excitations but also collec- 
tive excitations might be well-described with reasonable 
accuracy. At present, many theoretical techniques have 
been developed to solve mciidela.mi±h short-range interac- 
tions. The Hubbard modeOc^'ElP is a typical example, 
which is known to show various kinds of phase transi- 
tions and critical phenomena. Thus, if we can improve 
our technical scheme to unify DFT and knowledge of the 
strong correlation effects, our theoretical understanding 
of materials would be wider and richer. 
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We can say that the so-caUcd LDA+U (ancLipAj-ji^ 
approach has been developed along this line.BESEa'ESP 
In the LDA+U, the KS orbitals are used as a basis for an 
additional many-body problem. Indeed, idea of LDA++ 
is defined as a unified strategy in mhich the dynami- 
cal mean-field approximation (DMF)^^ is used to solve 
a many-body Hamiltonian, which is essentially an ex- 
tended Hubbard model given by the KS orbitalsES-* As 
well-known, DMF is a method to deal with short-range 
dynamical correlation correctly. However, the LDA+U 
energy functional is an artificial functional obtained by 
adding extra U-terms and subtracting a counter term, 
which is valid only in an atomic limit. Thus, it is desired 
to formulate a strategy to derive a proper energy func- 
tional by a mathematically reliable approach. As we will 
show, it is possible to determine the energy functional 
using a variational method which is same as that used in 
the Kohn-Sham theory. 

In this paper, we present a version of our formalism to 
obtain an extension of the Kohn-Sham equation. In this 
version, we consider a non-relativistic many-electron sys- 
tem. Thus, we will start from the energy functional given 
by Levyfl) We will introduce model functionals, which 
are composed of a kinetic energy functional, an interac- 
tion energy functional and a new exchange-correlation 
functional. A single particle density is given by a many- 
body state which is a solution of an auxiliary many-body 
problem. The energy functionals proposed in the present 
article are essentially extensions of the KS functionals. 

The original Kohn-Sham scheme has derived based on 
a next assumption. Namely, the true single-particle den- 
sity was assumed to be reproduced by a single Slater 
determinant and the Slater determinant was assumed to 
be a solution of an effective single-particle equation. This 
assumption is a so-called v-representability in KSE. To 
avoid usage of the v-representability, Hadjisavvas and 
Theophilou (HT) have introduced functionals on the 
space of Slater determinants instead of density func- 
tionals. Their strategy succeeded in bypassing the v- 
representability problem. However, their argument is a 
bit complicated for intuitive argument. Thus we will 
start our derivation using a variational method similar 
to the Kohn-Sham scheme from §^ A mathematically 
rigorous definition of our new functionals in accordance 
with the HT argument will be also presented in 

§2. Definition of the Problem 

We consider a non-uniform electron gas in an external 
scalar potential Wext(r)- Here the position vector r is 
defined in a three-dimensional space. The problem is 
given by the next non-relativistic Hamiltonian, 



V 



H = f + V+ /dV«ext(r)n(r) . 



(2.1) 



Using electron field operators with spin cr, which are de- 
noted as V'o-(r) and V'i(r), the kinetic energy operator T 
and the 2-body interaction V are given as follows. 



T 



2m 



d^r ^ lim 7/;^ (r')ArV'CT(i") 



(2.2) 



d^rd^r' 



7^t(r)V;t,(r')Va'(r')V'.(r) • 
(2.3) 

The operator, n(r) = X^o- V-'J(r)'0cr(r), gives the electron 
density, n{r), at r. 

The many-body problem given by H is difficult to solve 
in general. But we know that the system is stable, i.e. 
the energy is bounded below and that the ground state 
exists for any fext(r). Thus we introduce |<&gs) repre- 
senting the ground state of H. The single particle den- 
sity of I'I'gs) is given by. 



"■Gs(r) = (*Gs|fi(r)I3>Gs) 



(2.4) 



The rigorous version of DFT has been established by 
introduction of the next energy functional by Levy.EP 



F[n] 



min ($|r 



(2.5) 



Here, min$^„ represents the constraint minimiza- 
tion, where $ is optimized under a constraint that 
($|n(r)|<I>} = n{r). We know that the minimum for the 
Coulomb problem under discussion-.sxiats and that F[n] 
is a well-defined functional of ri(r)J3'QcP It is recognized 
that the strategy of the KS, approach should be based on 
the functional, eq. (p^.By# 



The central assumption to derive our extension of 
the KS scheme is the following. We assume that 
there exists an auxiliary many-body system and that 
the ground state of the model system, I^'gs), satisfies, 
(\E'Gs|'^(r)|5'Gs) — '^Gs(r). Our purpose is to find a sys- 
tematic strategy to determine the model system. This 
model should be described by a set of auxiliary equa- 
tions, which we introduce now. 

The first set of equations is to define single-particle 
orbitals. 



2m 



Ar + Woff(r) > 01, ^(r) = o-(r) 



(2.6) 



Here Vcsi'r) is an effective potential, which will be de- 
termined by a variational method. Spin-dependence of 
Woff(r) does not appear in our non-relativistic version, 
since we consider only the density functional and do not 
deal with a spin-density functional in the present paper. 
However, as we discuss in the last section, I^'gs) can be 
magnetic. 

We assume existence of Vcs{r) in the same man- 
ner as Kohn and Sham did. Here, we should note 
that the existence of Woff(r) assumed in the original KS 
scheme has not been proved nor disproved, although 
the v-representability assimied in the Hohenberg-Kohn 
theory was disproved. eIBcZI'' In we will utilize the 
HT technique to redefine our scheme avoiding the v- 
representabili ty p roblem. 

Since eq. ( p.q ) determines a complete set of single- 
particle wave functions, we can introduce Fermion op- 
erators, cl ^ and Ci^a- They create (or annihilate) the 
i-th state, 0i.(j(r), with a spin a and satisfy canoni- 
cal commutation relations for the Fermion operators, 
{c-i,<y T c^j fji} = Sijda^cr'. Note that we may use another 
basis set comprised of localized orbitals e.g. the Wan- 
nier basis, to consider an auxiliary many-body prob- 
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lem. However, the kinetic energy term becomes a tight- 
binding Hamihonian with hopping terms between every 
pair of orbitals. Thus, we utihze the eigen states of eq. 
(2.6) to show the formal strategy. We will soon explain 
the reason why the formulation is easily generalized using 
localized electron picture. 

Using cl ^ and Ci^^, the effective Hamiltonian of the 
auxiliary many-body problem is written as follows. 



(2) 



i,j ,k,l,(T,a' 

+ i?rxc . (2.7) 

Thus the second equation of the auxiliary problem is, 

iJeffI*) . (2.8) 

Again, effective interaction parameters, V^^\i, j, a) and 

(2) 

V^g {i,j,k,l;a,a'), should be determined by a varia- 
tional method. The following sections is devoted to ex- 
planation of a method to determine these effective in- 
teractions and the effective potentials in a sel f-consistent 
manner. The last term, i?ixc, in eq. (2.7) represents 
residual exchange-correlation interactions not described 
by two-particle interactions, V^g {i,j,k,l]a,a'). Defi- 
nition of Hj-xc will be given also using the variational 
method. 

An important point of our method is that a summation 
denoted by ^ is arbitrary as far as H^s is an Hermi- 
tian operator. Namely, we have degrees of freedom to 
optimize the auxiliary effective problem. In other words, 
we may consider Hubbard-type localized interactions at 
any time by considering a proper unitary transformation 
from extended states, (f>i^a{^), to localized orbitals. 

The many-body state, |^), should be given as a sum- 
mation of Slater determinants. We use a term, "a multi 
Slater determinant" as a summation of Slater determi- 
nants below. In the present notation, it reads, 

J 



/ = {il,(Ti,---,ljv,(Tjv} 



JO), (2.9) 



for an iV-particle system, where I represents a multi- 
index, {ii, CTi, • • • , lAT, ctat}, a nd c oefficients /; should be 
determined by solving eq. (2.8). Namely, we assume 



that the ground state of the second auxiliary problem, 
|^'gs)j is obtainable, when we derive the formal theory. 
We introduce the reduced density matrices. 

p{t,j,a)^mcl^cj.,,\^) , (2.10) 

p^^\i,3,k,l-a,a') = (*|ct^c],,,Cfe,,,Q,,|*) ,(2.11) 

{i, j, k, I; cr, cr') = p(^) k,l;a, a') 

- p{i,l,a)p{j,k,a') . (2.12) 

The single particle density of \^) is given as, 

n(r) ('l'ln(r)l'I') 



a 

(vE'|^c/,,*,(r)</.,-.(r)4,c,-.|vI/) 
E '^»>(i")0j,<T(r)p(i,j» . 



(2.13) 



Following the central assumption, 7i(r) — nGs(r), if the 
auxiliary problem is properly chosen and if l^f) — I^'gs)- 
The total energy of the auxiliary model is given as, 



(2) 



E V^l\i,j,k,l-a,cj')p^^\i,j,k,l-a,a') 



(2.14) 



§3. Model Energy Functional 

To find an extension of KSE, the second important 
step is determination of energy functionals for the auxil- 
iary model. We define a kinetic energy functional To[n] 
and an interaction energy functional V^el?^] as follows. 



To[n] 



(r)c: 



I*) 



2m 



0j,^(r)p(i, j, cr) 



E / '^^^<^»>W"<=ffW</'j,'T(r)p(i, j, cr) 
^eip(i,i,cr) - Y J d^rvcs{r)n{r) . 



(3.1) 



3' 



(2) 



d^r Sr' 



-n(r)n(r') 



i,j,k,l,<T,(T' 

X(t>k,a' (r')0/,<T (r)p^^^ (i, j, k, I; a, cr') 
1 



= E l^{^i.'y(t>j,a'\ — \<j)k,cr'(t'i,a)pii,l,cr)p{j,k,a) 



(2) 



i,j ,k,l,(T,a' 
(2) 



1 



k,a<l>j,cr'\ — \<Pk,a'4'l,<T)p^'^\hj, k, I; CT, Cr') 
k,a<l>j,a'\ — \(t)k,a'(l)l,cr){'^\cl„cl^^,Ck,a'Cl,a\'^) 
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(2)' 1 g2 

+ 2^'l>^,'^'^J,'^'\ — \'l>k,a'(l}l,a){'f\claCl,a\'f) 

i,j,k,l,(T.(T' 

x{^\cl^,Ck,.,m . (3.2) 
Here an interaction parameter is defined as, 



i'j.a'\ — \(t'k,cr' 4>l,a) 

r 



d\d?r' (r)0*,,(r')0.,.'(r')0i,.(r) , 

|r — r I ■' 

and a summation represents =y^ — y^}'^\ 

We have to caution that in both eq. ^'^^ 



(|3.2| ) I*) is imphcitly assumed to be I^'gs)- Namely, \^) 
is optimized with a condition, (\['|n(r)|4') = n(r). Then, 
the reason why T^ln] and K,o[n] are functionals of the 
density is explained as follows. Since orbitals, (/>i^cr(r), 
are decided by a KS-like equation, eq. ( ^.6[ ), these or- 
bitals are a functional of n(r). Since the second equation, 
eq.(p^), is determined by (j)i^„{v), |^gs) and resulting 
density matrices are also a functional of n(r). Details of 
the simulation process are also discussed in § |5.2| . How- 
ever, the above reasoning, which accords with the orig- 
inal KS argument, is not accurate enough. Namely, we 
have to specify a minimization procedure in detail. As 
clarified in this process can be regarded as a search 
in the space of n{v). Thus, let us utilize notations To[n] 
and Vi3e[?T.]. 

Another energy functional which we call the residual 
exchange-correlation energy functional, i<^rxc['^], is de- 
fined as the difference between the true energy functional 
and To[n] + V^cM as. 



F[n]^To[n] + V,M+E,^M 



(3.3) 



Thus the total energy functional for an electron system 
in an external potential Woff (r) is given as follows. 



E[n]=TQ[7i\+V,, 



- Er- 



+ / d^rweff(r)n(r) . (3.4) 



Here we should comment the next point. In the last 
expression of eq. (|3.2|), selected interaction processes are 
explicitly kept as 2-body interactions in X)^^^- But, for 
other processes, a corresponding Hartree term is used 
and corrections from remaining exchange-correlation ef- 
fects are included in _Bi.xc [n] ■ 

§4. Variational Method 

To determine effective interactions and effective poten- 
tials, we utilize a variational method. Consider variation 
of £J[n] with respect to variation of n(r). The variation, 
(5n(r), is separated into two parts. 

5n{r) = 5n'^°\r) + 5n'^'^\v) , (4.1) 

<5n(°)(r) = ^ {54>l,{v)<t>,A^) + 0*,(r)<50,,,(r)) 

xp(i,j», (4.2) 
<5nW(r) = ^lAr)<PjAr)Spit,J,a) . (4.3) 



These equations mean that the density variation can be 
taken with respect to two types of independent direc- 
tions in the functional space: one is by S(j)i^a{^) and 
the other is by Sp(i,j,a). To be more precise, we note 
that our model problem is solved by performing two suc- 
cessive steps alternately. The first one is to determine 
the single-particle orbitals given by eq. (|]^) fixing all 
density matrices and/or coefficients /; of the many-body 
state I^'gs)- After redefining (j)i^^{r), the many-body 
problem given by eq. (^|^) is solved and density matri- 
ces are rebuilt. Thus, every variation of each physical 
quantity is taken with respect to S6i^[r) and 6 p(i, j, a) 



Our system of equations, i.e. eq. (2.6) and eq. (2.8), has 
to be solved self-consistently as a whole, so that Sn^°\r) 
and Sn^'^^r) are mutually affected with each other. How- 
ever, to derive determining equations we can consider the 
first order variation for each quantity. Therefor, we can 
classify variation into two types, i.e. those related to 
orbital variation and others related to variation of the 
density-matrices which are denoted by (o) and (d). 



4.1 Variation of the Energy Functional with respect to 
the Single- Particle Orbitals 
In this subsection, we consider Sn^°\r) which is di- 
rectly given by 5(j)i^^{r). Before presenting the result, we 
derive some relations. Let us consider variation of eq. 
(2.6), where and Weff(r) are also dependent variables 
of Sn{r). 



-T^A + v^sir) )■ S(f),,^{r) + Jueff (r)(/>,;.<^(r) 
2m J 



(4.4) 



Multiplying eq. ( |4.4| ) by 4>l„{r) (or 0*^(r)) from the 
left, integrating the result with respect to r, we obtain 
the following useful identities. 



d\6vM\<l>,,Ar)\' 



(4.5) 



(£j- -£,) / dV0*^(r)J0,,^(r) 



= - / dV(5wcff(r)0*^(r)<?!),,<,(r) . 



(4.6) 



Variation of the total energy functional due to Sn'-°\r) 
is obtained as follows. 



+ y"dVz;e,t(r)<Sn(°)(r) , 
where StI^"^ [n] and SVc^A [ n] are given as, 

= ^(5ejp(i,z,cr) - Jd^r 5ueff(r)n(r) 
^dVueff(r)(5n(°'(r) 



(4.7) 
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'i,j,/c,/,cr,(7' 



/ L.J ,t\,,L.U ,U 
d^rcj)* {r)S(t>j^a{r)p{i,j,(T) f ( f 



dV«eff(r)(Sn(°)(r) , (4.8) --yeffW} (Jn^"*' (r) . (4.13) 

Thus the total variation of the energy functional is 
rTr(o)r i given as follows. 

(2) w 2 ^^[" 



^ 1 f 

= l^\{64>^,a4>3,a'\-\4>k,a-<t>l,a) = Y e,S p{l , i, a) 

i,(T 



+ (0j,o-'50j.<t'| — I'^fe.o-"/';,^-) + {4>i.a4'j.a'\ — \Hk.a'4ll,a) 
+ {(t>i,a4>3,a'\ — \(t>k,a'54>l,a) \ p'^^Hi , j , k , l] a, a') 



r 

o2 



(2) 



If 

+ / d^rd^r' " n{r')6rS°\v) . (4.9) ^.3M,L<y,'r' 



i'i,a<t)],a'\^\<l>k,a' <l>La)S p^'^^i , j , k, I; CT, Cr') 



Thus we obtain the next expression for SE'^°\ 



= Jd^r<j)*„{r)6<j).j,a{r)pii,j,<T) 



^^'^|y^V^n(r') + ^+.e.(r) 
. -«cff(r)}(5n(r). (4.14) 

(2) 2 This equation tells us that the energy functional becomes 

^ -S{4>i^c,<t>],a' I — l0fc,<T' '/';,ct) (i, j, fc, cr, a') stationary with respect to the density variation, if a next 
i,j,k,i,(7,iT' ^ ^ equation is satisfied for the solution of the auxiliary prob- 

/dV /dV^-n(r') + f^+z;ext(r) \- ^ 

J [J |r-r'| (5n(r) 2^erdp{i,i,a) 

-v,s{r)} 6n^°\r) . (4.10) 

('^^ "'^j) / d^r(j)l„{r)6(j)j,air)pii,j,(T) 



4-2 Variation of the Energy Functional with respect to ^ 
the Density Matrices 
The variation Sn'^'^\r) is given by variation of the den- 
sity matrices, Sp{i,j,a) and Sp^^^{i,j,k,l;a,a'). Varia- 
tion of energy functionals, (5rQ''''[n] and (5V^lf'' [n] , due to ' ' 



k.l,cr.(T' 

<5nW(r) are given as, " _^ (^...^...H- J, fc, a')| = (4.15) 



r 

As shown in the next subsections, this condition for the 



= ^ei(5p(i,i,cr) - / dVwoff(r)(5n('')(r) , (4.11) density matrices is satisfied, if |*) is a solution of a many 
■J bodv nroblem sfiven bv en. (U^) with nronerlv choser 



body problem given by eq. (^_^) with properly chosen 
effective interactions. In other words, the present sta- 
tionary condition on E[n] is guaranteed provided that 
we can determine effective interactions such that the so- 



2. lution of eq. ( |2.8| ) satisfies eq. (4.15). In the foUow- 

) ing subsections, we show existence of required V(.g{r), 



M.M.-.-' - V}^\ij,a), V}^\i,j,k,l;a,a') and H, 



e2 



+ J d^rd^r' y——^n{Y')5n^'^\Y) . (4.12) ^ g Variation of the Second Auxiliary Equation 

Now we investigate variation of eq. (2.8). We note 
We obtain the next expression for 5E'^'^\ that not only the solution |*) but also aU parameters, 

5E^'^\n\ ^^"^ orbital energy, interaction parameters and -ffrxc. 
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are varied. 

i.(T 
(2) 

+ (i?,,c+^i?rxc)}|* + <5'J') 

= (£; + (5£;)|«' + ^^') . 



(4.16) 



As a result, we have the next equation which holds in 
the first order of variation. 

5E 

= '^{6eip{i, i, cr) + e^Spii, i, a)) 

',fT 

(2) 

E I'^^off^ ihJ, k, I; a, a')p('^ (z, J, fc, I; a, a') 



+ V^l\i,3. k, I; a, a')Sp^^Hi,j, k, I; a, a')} 



(2) 



(2) 



+ J d^rSvcs{r)n{r) 

(2) 

+ E {^^csi'^^j^k,l;a,a')p{i,l,a)p{j,k,a') 

i, j,k, I, <t,(t' 

+V^g'{i,j, k, I; cr, cr') cr)p(j, /c, cr') 

+p(i,/,cr)(5p(j,A:,cr'))} 

+ E {^Kff^(^,J,^)p(*,J,^) + Kff^(^,J,^)M*,J,^)} 



(4.17) 



Here we have utilized eq. ( |4.5|) and eq. (4^). For the 
stationary solution of eq. ( p.§| ), SE = and we have 
another useful relation. 



4-4 Determination of Effective Interactions 

The expression of eq. ( 4.14 ) tells that natural expres- 
sion for VcS (r) is similar to the KS effective potential and 
reads as, 

d'r ^^n(r') + "—^ + i,ext(r) . (4.18) 
|r-r'| dn[r) 

Then, we have an expression of (5woff(r) as, 



Svcs{r) = / (fr -—-—-Sn{r') + 



6n(r')Sn{r) 



Sn{r') 



(4.19) 

Besides, comparing eq. ( 4.14 ) and eq. ( 4.17 ), we notice 
that V^g {i,j, k, I; cr, cr') should be chosen as. 



1 e 

V}pii,j,k,l;a,a') = -(0.,.0j,.' I^l'^fc,^'^/,^) • (4-20) 



Substitute eqs. ( [4.20D into eq. ( [4.171 ) and let SE = 0, 
then we find that the required stationary condition, eq. 
( 4.15| ), is satisfied, if the next condition holds. 



E {sv}^\^,J,^)pi^,J,<^)+v}^\t,j,a)6pii,j,a)] 



-(^(*|i/„c|*) 

tPr (5wcff (r)n(r) 

(2) 



- E [svlsihj^k,l;a,a')p{i,l,a)p{j,k,a') 

i,j ,k,l.(T,a' 

+Ks («' k, I] cr, a') {Sp{i, I, (j)p{j, k, cr') 

+p{iJ,a)Sp{j,k,a'))} (4.21) 

This new condition is useful for determination of 
V}ff\ i,j,r) and_Hrxc. Substitute eq. ( 4. IE ) and eq. 
(4.20|) into eq. (4.21), we obtain the next condition for 



V]g'{i,j,a) and iJ^xc- 

E {Sv}^\i,j,a)p{i,j,a) + V}^\i,j,a)6p{i,j,a)} 



+ 5(*|i?,xc|*) 



x(50*,(r')0,-4r') + 0*,(r'),50,-,(r'))p(*,j, a) 
E / d'^'^'^' ^^«(r)0*.(r')</>,..(r')<5p(z,j,a) 



(2) ( (2) 



i 1 e 
E 2^i<Pt.'r(t^k,a'\ — \(l)l,a'(f>j,a)p{k, I, <j')p{i,j, u) 
k,l,a' 

E i4>i.<y't'k,c' I — \(t)Lcr'4'j,<y)pik, I, <T')Sp{i, j, Cr) 



d^rd^r' 



-n{r)6n{r') 



(4.22) 



(5n(r')(5n(r) 

Here, a summation, represents a restricted sum 
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mation satisfying YlUa T,kla' = T,ik,i,j,a,a' ■ The con- 
dition, eq. ( 4.22| ), is satisfied, if we choose V^ff\i,j,r) 
and HrTcr as follows. 



d rd 



(2) 

(2 



X > ^[(Pi,a(f>k,a'\ — \(l)l,a'(l)j,a)p{k,l,a'), (4.23) 



SE 

(*|i/rxc|«') = - jd?r -^^n{v) + S,,e W , (4.24) 



(5ri(r')(5n(r) 



n(r)(5n(r') 



(4.25) 



Here we note that required conditions for //rxc are only 



eq. (4.24) and eq. (4.25). These equations are satisfied, 
if we choose an explicit expression of -ffrxc a-s, 



(4.26) 



SE 

ErxcN = - / ^^ri(r) + SrxcM . (4.27) 

The operator, /, represents an identity operator in the 
phase space of the A'^-particle system under discussion. 
Finding of these expressions for effective interactions is 
an important result of the present theory. 

§5. The System of Auxiliary Equations 

We have obtained a new set of equations, which gives 
a generalization of KSE. The system is summarized as 
follows. 



-T^Ar + Weff(r) > <P^,a{r) = e,(j),^^{r) , (5.1) 
2m ' 



(2) 



}\^)=E\^) , 



-off (r) - J d\ ^«(r') + ^ + «cxt(r) , (5.3) 



V^l\i,j,k,l\a,a') = ]^{(t)t,a4>j,a'\^\4)k,a'4>i,<y) , (5-4) 



(2) 



ll 1 



k.l.a' 



ErxcM = £^rxcM " /rf^?" ^^^^(r) . (5.6) 

J dn[v) 

In the following subsections, we discuss interpretation 
of the new equations and possible calculational strate- 
gies. 

5. 1 Interpretation of the Equations 

Each equation in the new set of auxiliary equations 
has its own meaning. We review important properties. 

Eq. (5.1) is a direct generalization of KSE. The effec- 
tive single-particle potential, eq. (5.3), is comprised of a 
Hartree term, a new exchange-correlation potential and 
the external potential. However, we have to comment 
that n(r) is given by a multi Slater determinant, |^'gs)j 
and definition of E'rxc is different from that utilized in 
the original KS scheme. 

Eq. (5.2) gives an auxiliary many-body problem. The 
summation, is arbitrary as far as HcS is Her- 

mite. If X)*'^^ is taken for all possible combinations, 
E'lxc = and the problem becomes identical to the orig- 
inal Coulomb-potential problem for a non-uniform elec- 
tron gas. If the summation is not taken for any combi- 
nation, our system of equations reduces to KSE. Thus, 
our new equations is a natural extension of KSE which 
makes a connection between the usual KS scheme and 
the Coulomb problem. 



We have obt ained eq. (5^) from observation of eq. 
( 4.14 ) and eq. (|4.17 ). A natural question may be "Why 
does any screened interaction appear?" This point is 
discussed in §^ too. 

Eq. ( |5.5| ) is naturally interpreted as a counter term 
required by definitions of eq. (5^) and eq. (5^). This 
term represents a Hartree-like mean-field. However, we 
have to be careful with appearance of the second term 
with a sum mation X) • 

Eq. (|]^) is a correction coming from the reduced ex- 
change correlation functional. This expression indicates 
that Erxc is a functional of n(r). 



(5-2) 5,2 Calculational Strategy 



If we assume that we can solve eq. ( p.lD and eq. (^) , 
a process to obtain the ground state will be the following. 
Here, we also assume that we know a value of -Brxc [n] and 
its variational derivative for any n(r). 

1. Start from an initial single-particle density n[v). 

2. Once n(r) is gi ven, 'yoff(r) is given. 

3. Solve eq. (|5.l| ) to obtain (/)i.o-(r). Here ri(r) is as- 
sumed to be fixed. 

4. Solve eq. ( ^ ) by a proper method. Since 
^as' J' -^rxc are given by the density matri- 
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ces and n(r) of the solution, the determining equa- 
tion has to be solved in a self-consistent manner. 
5. Recalculate n(r) and goto the step 2 until conver- 
gence in n{r) is obtained. 

The above optimization process is done in a phase space 

of n{r) given by 

§6. Interpretation of New Energy Functionals 

As discussed in we considered that functionals, 
ro[?^]i K!c[?^] and Si-xcM, are density functionals. We 
discuss this important point in more details. We hope 
that the explanation below will answer remaining ques- 
tions for our formalism of DFT. 

Let us review the HT theory at first. To formulate 
rigorous scheme of KSE, they considered a functional of 
a single Slater determinant Here we use a notation, 
n^r) = (0|n(r)|0). 

Gh[<I>] = {cl)\f\cl)) + AT{<j>) + Jd\v,^t{r)n^{r) 
If 

+ - Jd^rd^r' ^——^n4r)n4r') 

+ E^,{4>) . (6.1) 

To define AT((/)), we are required to determine a multi 
Slater determinant IV-") and a single Slater determinant 
\(j)') via two minimization processes given by next two 
functionals. 



Gt+v[<P] = min {i;\T + V\^) , 
GtW\ = min (0'|f 10') . 
Using IV'} and 10'), AT is defined as, 

Ar[0] = (v|f|v)-(0'|f|0'). 

Another functional -Exc(0) is given as, 



(6.2) 
(6.3) 

(6.4) 



i?xc[0] = (V-lv^lV-) 



If 

2 / irr^"'^^'")"'^^'"') 



(6.5) 



Since Gh is rewritten as. 

Gum = imc^) - ('^'iri0') 

+ {^P\f + V\^P)+[ dVwext(r)n4r) , (6.6) 



it is easily shown that the minimization of Gh[4'] with 
respect to |0) is equivalent to obtain the true ground- 
state energy, Eq. Namely, 



min Gh\4'] — min < F[n] + 



^^;ext(r)f^(I■) j ^ Eq , 
(6.7) 

where F[n] is Levy's DFT. The minimization of Gh 
respect to |0) is shown to be identical to solve KSE.E3' 

The strategy given in the present paper is to utilize a 
multi Slater determinant |\1/) in place of \(f>) and replace 
T with T-t-K-educed defined below. Namely, our extension 
of KSE is rigorously given by the following functional of 
1^*). We use a notation, n,p(r) = ('I'|n|^'), too. 



= (*|T + K-cd|*) - min (*'|T + K-od|*') 
+ F[n^] + / d^rt;oxt(r)n*(r) . 



(6.8) 



Here a reduced interaction operator, V^od, is defined as a 
following operation. 

(^iK-cdl^-) 

(2 



r — r' 



- (vI/|4,Q,.-|*>(*|c,Vcfc..,|*)} ■ (6.9) 
Then, we can rewrite Gh as. 



Gnm 



(^-iT+K-edl*) 



d^rwcxt(r)n* (r) 



+ / d^rd^r' ^ n^(r)n^(r') 



r r' 



+Er^c[n^i,] ■ 
Here -El■xc[^^*] is defined as. 



(6.10) 



= min (^'iT-fV^edl*') 

-If d\d\' ^-^n^{r)n^{r') . (6.11) 
2./ r - r 



This definition of -Erxc is equivalent to eq. ( |3.3D , if we 
identify nq,(r) with n(r). It is straight-forward to derive 
our set of auxiliary equations by a variational principle 
on functionals Gh- 

These expressions tell a lot about possibility of our 
extensions. Eq. ( |6.9| ) requires to determine not only 
range of X)'^' t>ut also KS orbitals, (/)i^o-(r). However, 
we may utilize any localized orbitals in the definition of 
Vied- Since eq. ( |6.9[ ) takes a form of a correction to 
the mean-field approximation, we should include a key 
process, e.g. interactions between electrons in the same 
localized orbital, to reduce amount of E^xc- Besides, we 
might be allowed to utilize a screened interaction in eq. 
(6.9), if -Erxc became a small correction. 

Note that definition of Vred requires to fix interaction 
terms defined by ^^^^(r). In the previous derivation, we 
implicitly assumed an optimization process of orbitals, 
0i^CT(i"), too. Once 0i,cr(r) is varied, values of E^^c change 
in principle. In several expected situation like the valence 
fluctuation, however, the Hartree-term would be dom- 
inant to determine the shape of local orbitals. Thus, 
to approximate -Erxc independent of a small change of 
0i^cr(r) might be practical. 



Gnm 
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§7. Possible Approximation 

For a practical calculation, we have to decide an ap- 
proximate evaluation of i?rxc as the original KS scheme. 
However, we have to consider that -Erxc depends on the 
summation, in eq. (5^). As a case study, let 

us consider a localized f-orbital embedded in conduction 
electrons. In this discussion, we omit the relativistic ef- 
fects for simplicity. 

Using a proper unitary transformation, a Hubbard- 
type local interaction between f-electrons can be de- 
scribed explicitly by 2-body interaction terms in eq. 
(5.2). Let us assume that other 2-body interactions 
are excluded in ' . This effective Hamiltonian corre- 
sponds to an Anderson impurity HamiltonianJ^S-* where 
both the f-orbit al a nd conduction states are given by so- 
lutions of eq. (5J). The residual exchange-correlation 



term, i?rxc, represents not only all of the exchange- 
correlation effects among conduction electrons but also 
a part of exchange-correlation effects between an f- 
electron and conduction electrons. The latter exchange- 
correlation includes an effect by charge fluctuation on the 
f-level, which introduce effective expansion or contrac- 
tion of f-charge affecting the effective potential, Voff(r), 
too. However, we can expect that most of the effects 
coming from this charge redistribution can be repre- 
sented by the Hartree term in Vcs{r). 

Thus, we can expect that -Erxc is not large in amount, 
as far as we include the most important exchange- 
correlation processes in our auxiliary equations through 
the local interaction, V^g {i,j,k,l;a,a'). The simplest 
practical approximation for i?rxc is a kind of the local 
density approximation (LDA). Many strategy of LDA 
for -Brxc would be possible. However, to determine this 
function, we have to decide inputs. If we adopt repre- 
sentation by localized orbitals,E3' we can define a single 
particle density of f-electrons, fT.f(r), and that of con- 
duction electrons, nc(r). Then, some versions of LDA 
would be classified by the inputs as, 1) n{r) at r, 2) rif (r) 
and ric(r) at r, 3) nf(r), nc{r), Arif(r) and Anc(r) at r. 
The categoot^S) corresponds to GGA or more closely to 
meta-GGA.EH' Determination of a practical method is a 
remaining important problem. 

Reliable ii^jxc should be determined by accurate nu- 
merical calculations. As has been done for determina- 
tion ^jf_LDA using the diffusion Monte-Carlo (DMC) 
dataMEiP for the uniform electron gas,EJ parameteriza- 
tion might be possible if there were data for an impurity 
problem by DMC or other accurate methods. 

Another possibility is to adopt the optimized effec- 
tive potential method witja-esact evaluation of the ex- 
change interaction (EXX)^'^^'^^ For example, if eval- 
uation of i?rxc byi— pXX with RPA for correlation en- 
ergy (EXX-|-RPA)Ea) is tractable for evaluation of iJixc, 
and if the short-range dynamical correlation described 
in Hcs is well-described by DMF, a hybrid method of 
DMF+EXX+RPA might be possible. 

§8. Discussion 

Starting from Levy's DFT, we have formulate an ex- 
tension of the Kohn-Sham equations in order to simulate 



correlation effects coming from localized nature of elec- 
trons in materials. The system of equations consists of 

1) determination equations of single-particle orbitals and 

2) a many-body model problem. These auxiliary equa- 
tions should be determined self-consistently. In §^ we 
have given a rigorous derivation of the extended KSE. 

As commented in § ^.1| , our new scheme provides a con- 
nection of the original Coulomb problem and KSE. This 
property is very useful for application. Conversely, a se- 
rious problem is how to determine most efficiently. 
Without ruining qualitative accuracy, we have to control 
this summation to optimize simplicity for actual calcula- 
tions. The problem, we think, may be solved by knowl- 
edge of information science. 

From a practical point of view, our theory may lead 
to a simple approximation method including non-local 
correlation effects, which cannot be treated even by 
GGA_l3' Comparison between our formalism and meta- 
GGAE3'' might be an interesting problem. Since our 
strategy supplies a kind of the so-called first-principles 
calculation, we hope that the present scheme may be a 
step toward a unified theory for the interacting electron 
systems. 



Eq. (5^) can be interpreted as a kind of the extended 
Hubbard model. Although 2-body interactions appear- 
ing in eq. (5.2) is bare Coulomb (and exchange) inter- 
actions, the model system can describe every magnetic 
systems as well as non-magnetic materials. Namely, if we 
properly solve the auxiliary problem, spontaneous sym- 
metry breaking can be discussed in a realistic electronic- 
state calculation. Besides, if a superconducting state was 
realized only by electron-electron repulsion, our strategy 
could be utilized to explore the superconductivity by a 
DFT-based calculation. |_. |_. |_. |_. |_. 

As well-known, the relativistic DFT (RDFT) NN piMkI) 
the current density functional theory (CDFT)c3'c3 and 
the unified version, i.e. the_rclativistic current density 
functional theory (RCDFT)c3'La) have been established. 
There should hopefully appear an extended version of 
the present theory based on RCDFT. Besides, there are 
many remaining important directions to be investigated, 
e.g. an extension of KSE for correlated systems at finite- 
temperaturesEK_and_DFT-description of time-dependent 
phenomenaEO' o cZfCa^ 
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